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Abstract 

We define the local periodic linking number, LK, between two oriented 
closed or open chains in a system with three-dimensional periodic bound- 
ary conditions. The properties of LK indicate that it is an appropriate 
measure of entanglement between a collection of chains in a periodic sys- 
tem. Using this measure of linking to assess the extent of entanglement 
in a polymer melt we study the effect of CReTA algorithm on the en- 
tanglement of polyethylene chains. Our numerical results show that the 
statistics of the local periodic linking number observed for polymer melts 
before and after the application of CReTA are the same. 



1 Introduction 

Polymer melts are dense systems of macromolecules. In such dense systems the 
conformational freedom and motion of each chain is significantly affected by 
entanglement with other chains which generates obstacles of topological origin 
to its movement [1] [2] . These obstacles are called topological constraints (TCs) 
(Figure P). Edwards suggested that the local interaction between the chains 
restricts their motion to a tube-like region. The axis of the tube is a coarse- 
grained representation of the chain and it is called the primitive path [I] . Owing 
to its simplicity, the tube model allows us to carry out a theoretical analysis and 
describe the dynamical and rheological properties of polymer melts [3 El HI EE]. 
The CReTA (Contour Reduction Topological Analysis) algorithm [7] employs 
the Doi-Edwards perspective to study the entanglement in polymer melts. It is 
a coarse graining algorithm that reduces a computer generated atomistic sample 
to a simpler entanglement network of primitive paths (Figure [lj , where inside 
the vertices occur TCs. 

The computer simulation of a polymer melt usually involves molecular dy- 
namics and Monte Carlo simulations [H |H1 QIIj ■ In order to eliminate boundary 
effects, for the modeling of a bulk system Periodic Boundary Conditions (PBC) 
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Figure 1: (a) Representative atomistic PE sample and (b) the correspond- 
ing reduced network, (c) The central cell C and the periodic system it gen- 
erates. The generating chain i (resp. j) is composed by the blue (resp. 
red) arcs in C. The free chain / (resp. J) is the set of blue (resp. red) 
chains in the periodic system. Highlighted are the minimal unfoldings of the 
images I\ and J\. For the computation of LK{J,I) we have LK(J,I) — 
L(Juh) + L(J\,h) + L(Ji,I 3 ) + L(Ji,h). We can see that L(I\,J\) = 
LOA.iiJ.LOTi.Ja) = L(Ji,/ 4 ),L(/i, J a ) = L(J U I 2 ) and L{h,J A ) = L(J U I 3 ), 
thus LK(I, J) = LK{J,I). 



are applied [11] , In a PBC model, the cubical simulation box or cell is replicated 
throughout space to form an infinite lattice. In the course of the simulation, 
when a molecule moves in the simulation box, its periodic image in each one of 
the replication boxes moves in exactly the same way. Thus, as a molecule leaves 
the simulation box, one of its images will enter through the opposite face at ex- 
actly the same point. The configuration of the molecules inside the simulation 
box generates the entanglement present in the continuum. Because of this, the 
periodic system induces special features in the entanglement of the simulated 
polymer melt. For each chain in the melt we can distinguish between two types 
of topological constraints imposed on it by the other chains: local topologi- 
cal constraints, which can be observed inside a simulation box, and large scale 
topological constraints, which are observed only at a larger length. 

The purpose of this paper is to describe methods by which one may quantify 
and extract local entanglement information in a system with PBC. More pre- 
cisely, in Section 2 we describe the structure of the PBC polymer melt model. 
In Section 3 we introduce the periodic linking number and the local periodic 
linking number as measures of entanglement in such a system with PBC and we 
study their properties. In Section 4 we apply the local periodic linking number 
to polyethylene (PE) melts and study the effect of the CReTA algorithm on the 
linking present in the system. 
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2 Periodic Boundary Conditions 



We study a polymer melt consisting of a collection of polymer chains of length 
to, by dividing space into a family of cubic boxes of volume L 3 , where L is the 
length of an edge of the cube, so that the structure of the melt in each cube 
is identical, i.e. we impose PBC on the melt Specifically, we make the 
following definition : 

Definition 2.1. A cell consists of a cube with n arcs embedded into it such that 
arcs may terminate only in the interior of the cube or on a face, but not on an 
edge or corner, and those arcs which meet a face satisfy the PBC requirement. 
That is, to each ending point corresponds a starting point at exactly the same 
position on the opposite face of the cube. See Figure \V[b) for an illustrative 
example. 

A cell generates a periodic system in 3-space by tiling 3-space with the cubes 
so that they fill space and only intersect on their faces. This allows an arc in 
one cube to be continued across a face into an adjacent cube and so on. Since 
physical chains are of finite length, we require that resulting chains must be 
compact. For practical purposes, only a finite number of copies of the cell are 
considered. We call such a collection a finite periodic system. 

Without loss of generality, we choose a cell of the periodic system that we 
call generating cell. Then any other cell c in the periodic system is a translation 
of the generating cell by a vector c — (c x ,c y ,c z ),c x ,c y ,c z G L7L. A generating 
chain is the union of all the segments inside the cell the translations of which 
define a maximal connected arc in the periodic system. For each arc of a gen- 
erating chain we choose an orientation such that the translations of all the arcs 
would define an oriented arc in the periodic system. For each generating chain 
we choose without loss of generality an arc and a point on it to be its base 
point in the central cell. The smallest union of the copies of the cell needed 
for one complete unfolding of a generating chain shall be called the minimal 
unfolding. The smaller number of copies of the cell whose union contains the 
convex hull of the complete unfolding of a generating chain shall be called the 
minimal topological cell. 

For generating chains we shall use the symbols A generating chain 

is said to be closed when the connected component that its segments construct 
in the periodic system is a closed chain. The collection of all translations of 
the same generating chain i shall be called free chain, denoted /. A free chain 
is identified with the collection of its connected components, each of which is 
a translation of an unfolding of a generating chain in the periodic system. For 
free chains we will use the symbols I, J, ... . An image of a free chain is any 
connected component in it. For images of a free chain, say /, we will use the 
symbols 11,1%,.... A free chain is said to be closed, when its components are 
closed curves. 

Proposition 2.2. Given a generating chain i in a cell, and its corresponding 
free chain I in the periodic system, every image of I contains a translation of 
each arc of i exactly once. 

Sketch of Proof. Let i' m , i' m+ i, ■ ■ ■ , i' m+ i denote a sequence of arcs that we meet 
along I\. Each one of these arcs is a translation of an arc of i. Suppose that 
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Kn+i = Kn- Then 5t must be that C+i+i = C+l> g ivin § C+2J = C- Thus J i 
not compact, contradiction. 

□ 

Proposition 12.21 shows that each image contains exactly one base point. We 
call the image of / whose base point lies in the generating cell the parent image 
and any other image of / can be defined by a vector with respect to the base 
point of the parent image. 

The geometric and topological information resulting from a periodic system 
is carried in its generating cell, but the global entanglement structure seen in the 
periodic system may be different from that seen inside a cell. Local topological 
constraints (local TCs) are those tight conformations of the arcs inside a cell 
that form barriers to small movements of a chain since a chain is prohibited from 
passing through itself or any other chain inside the cell by motions keeping the 
endpoints of the arcs of the generating chains fixed. 

3 Linking in PBC systems 

The Gauss linking number is a traditional measure of the algebraic entanglement 
of two disjoint oriented circles that extends directly to disjoint oriented open 
chains 

Definition 3.1 (Gauss 1877). The Gauss linking number of two disjoint ori- 
ented curves l\ and I2, whose arc-length parametrization is 7i(t), 72(s) respec- 
tively, is defined as a double integral over l\ and 

4W [0 ,i] i[o,i] I71W-72OOI 3 

where (71 (i), 72(3), 7i(i) — 72(5)) is the triple product of 71 (i), 72(5) and 71 (i) — 
72 0). 

In the case of closed chains the Gauss linking number is a topological invari- 
ant. If it is equal to zero, the two chains are said to be algebraically unlinked. 
For open chains the Gauss linking number is a continuous function in the space 
of configurations and as the endpoints of the chains tend to coincide, it tends 
to the linking number of the resulting closed chains [13] . For open chains, the 
Gauss linking number may not be zero, even for chains whose topological cells 
do not intersect. But as the distance between them increases, the Gauss linking 
number tends to zero [2]. We note that when one considers open chains, the 
concept of topological linking does not apply, as open chains are always topo- 
logical^ unlinked, i.e. they may be deformed so as to lie in disjoint topological 
cells. 

In the case of a periodic system a different measure of entanglement is needed 
in order to capture all the TCs and the effect of the periodicity of the confor- 
mations. Here, we apply the Gauss linking number and its extension to open 
chains to the situation of chains defined in a PBC context providing a measure 
of the large scale entanglement between the chains. We propose the following 
definition of the linking number of two free chains in a system with PBC: 
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Figure 2: An example for which LK ^ LKp. The red curve does not intersect 
the minimal unfolding of the blue curve, but they are topologically linked. If the 
blue chain moves into the blue cell, then the red curve will intersect its minimal 
unfolding and LK = LKp. 



Definition 3.2 (Periodic linking number). Let I and J denote two free chains 
in a periodic system. Suppose that I\ is an image of the free chain / in the 
periodic system. The periodic linking number, LKp, between two free chains / 
and J is defined as: 

LK P (I,J)=Y,L(IuJu) (2) 

u 

where the sum is taken over all the images J u of the free chain J in the periodic 
system. 

For the purposes of applications to polymer melts, where the density is rel- 
atively high, the topological constraints have a strongly local character and one 
may wish to focus attention on the local TCs [T51 El H51 E] - As a consequence, the 
following local periodic linking measure may provide a more relevant measure 
of the entanglement observed in a polymer melt. 

Definition 3.3 (Local periodic linking number). Let I and J denote two free 
chains in the periodic system. Let Ji, J2, . . . , J n denote the images of J that 
intersect the minimal unfolding of an image of /, say I\. The local periodic 
linking number, LK, between two free chains / and J is defined as: 

LK{I,J)= L{h,J u ) (3) 

l<u<n 

Example 3: Consider the cell as shown in Figure[TJ:. Then for the computation 
of LK(I, J) we need to consider the images of J which intersect the minimal 
unfolding of I\, as shown in Figure [lb: LK(I,J) — L{Ji,J\) + L{Ji, J 2) + 
L(h,J 3 ) + L(h,J 4 ). 

We note that if the chains do not touch the faces of the cell, and thus lie 
entirely inside the cell, LKp and LK are equal to the Gauss linking number of 
the two chains. We shall focus on LK. In a sequel paper we will study LKp 
and its relation to LK. 

Remark 3.4. The local periodic linking number, LK, depends upon the size of 
the cell. Also, as the minimal unfolding may change under the motion of /, then 
LK may change as well. Thus LK is not a topological invariant for two closed 
free chains. The minimal unfolding of an image, say I\, is homeomorphic to a 



3-manifold, say M. If the genus of M is greater than 0, there may be an image 
of J that is linked with I\ but that does not intersect the minimal unfolding of 
I\ (see Figure [5]) . We note though, that the size of a simulation cube that is 
appropriate for the simulation of polymer melts and other polymer systems, is 
such that it is impossible for conformations such as the one shown in Figure [2] 
to arise. So, for all practical purposes, the local periodic linking number is an 
invariant under the isotopic motion of closed chains. 

3.1 Consequences and properties of the local periodic link- 
ing number 

Proposition 3.5. The local periodic linking number LK between two free chains 
I and J of a system with PBC is symmetric, that is LK(I, J) = LK(J, I). 

Proof. Let us consider an image of J, say J±, that intersects the minimal un- 
folding of an image of I, say I\. Then every other image of J that intersects the 
minimal unfolding of I± is a translation of J\ by a vector v^. Then we have that 
1\ intersects the minimal unfolding of each J u = J\ + u u . Thus the translation 
of 1 1 by — , i.e. I u = 1\ — v* u , intersects the minimal unfolding of J\. 

□ 

The local periodic linking number is a continuous function. Indeed, the set of 
points of discontinuity in the space of possible configurations is of measure zero, 
so LK is a continuous function of the chain coordinates almost everywhere. The 
following properties follow from the fact that LK is a measure of entanglement. 

(i) In the case of closed chains LK is an integer. Indeed, it is equal to half the 
algebraic number of crossings between a projection of an image of /, say 1%, and 
the projection of the images of J which intersect the minimal unfolding of I\ 
in any generic projection direction. For closed chains whose minimal unfolding 
has genus zero, then LK is invariant under movement of the chain. 

(ii) In the case of open chains, LK is equal to the average algebraic number of 
crossings between a projection of an image of /, say I\, and the projection of 
the images of J that intersect the minimal unfolding of I±. Moreover, as the 
endpoints of the generating chains move closer, LK tends to the local linking 
periodic number of the closed free chains. 

Further, the local periodic linking number has the following properties with 
respect to the structure of the cell: 

(i) LK captures the contributions of all the local TCs that the one free chain 
imposes on an image of the other, but it also takes into account the global 
geometrical structure of the chains. 

(ii) LK is independent of the choice of the cell and of the translated image of 
a free chain / that we use for the computation of LK(I, J) with another free 
chain J. 

(iii) For closed chains in a periodic system, one may amalgamate cells to create 
larger generating cells in a new PBC such that LK = LKp, and then LK is a 
topological invariant. 
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3.2 Advantages of the local periodic linking number for 
the study of the linking in PBC 

Systems which employ PBC consist of an infinite number of translations of 
a finite number of chains. Applying a traditional measure of entanglement 
would imply computations involving an infinite number or, at least, a very large 
number of chains. Furthermore, because of the periodicity of the system, there is 
a scale at which one has homogeneity, so that entanglement can be characterized 
by the local structure and the degree of its homogeneity. 

The local periodic linking number, LK, is computed using a finite number of 
cells and, therefore, a finite number of chains thereby reducing the complexity 
of its calculation. 

Ideally, one would like to compute a linking measure directly from one cell, 
but the arcs of the generating chains inside the cell are relatively short. In 
order to capture the greater degrees of entanglement, or even complex knotting, 
a large number of arcs must be employed in the creation of a complex chain. 

The local periodic linking number is computed between long polymer chains, 
taking into consideration the entire geometrical structure of the chains. More- 
over, if one wanted to constrain the analysis to a smaller number of cells, it would 
be impossible to give a definition of linking that is symmetric and captures all 
the entanglement present in a cell. Indeed: 

Proposition 3.6. The local periodic linking number LK(I, J) between two 
closed or open free chains I and J takes into consideration as few as possi- 
ble cells required to assess all the local TCs that the one free chain imposes on 
an image of the other and at the same time to be symmetric. 

Sketch of Proof. The minimal unfolding (mu) of I\ contains all the local TCs 
imposed to I\ by J. But we have that LK mu (I, J) = J2iL(h,ji) J2 n LiJi^n) 
= LK mu (J, I), where ji, i n are the translated arcs of j and i in the minimal un- 
folding of Ii and Ji respectively. By adding the missing terms in each summa- 
tion until we get a symmetric definition of linking we obtain LK mu = LK. □ 

4 Application to polymers 
4.1 The CReTA algorithm 

To study the entanglement complexity present in a given polymer system, one 
has to coarse grain the polymer chains at a level where certain geometrical char- 
acteristics relevant to entanglement become evident [17] ■ The CReTA (Contour 
Reduction Topological Analysis) [7j algorithm fixes chain ends in space, and by 
prohibiting chain crossing, it minimizes (shrinks) simultaneously the contour 
lengths of all chains, until they become sets of rectilinear strands coming to- 
gether at the nodal points (TCs) of a network of primitive paths. When there 
are no possible alignments left, then we shrink the diameter of the beads of each 
chain and continue the same process, until a minimum thickness is achieved and 
no alignment moves are possible. Figure [T] shows an atomistic polymer sample 
and the corresponding reduced network. 
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A LK B D=LK(org)-LK(red) 




C LK D D=LK(org)-LK(red) 

Figure 3: Figure A shows the normalized probability distribution of LK for the 
original and the corresponding reduced chains in a PE-1000 melt and Figure B 
the probability distribution of the corresponding differences. Figure C shows the 
normalized probability distribution of LK for the original and reduced chains 
after end-to-end closure and Figure D shows the corresponding differences 



4.2 Numerical Results 

The CReTA algorithm has been shown [7J [T3] to give a meaningful representa- 
tion of the underlying topology of a polymer melt by comparing it to experimen- 
tal data. Using the local periodic linking number LK , we determine the extent 
to which the CReTA method preserves the linking measure of entanglement in 
PE melts, or if critical entanglement information is lost in the reduction process. 
Also, we perform the same comparison after end-to-end closure of the chains of 
a PE melt and its CReTA reduced melt. 

As we mentioned earlier, linear polymer chains are not linked or knotted in 
the topological sense. A method to study the entanglement of linear polymers 
is to perform closure of the chains and compute the resulting knot type [181 119j . 
Using LK we will test if the method of direct end-to-end closure [3 [T51 l2"0"]is 
reliable and if LK can detect topological differences between the open and closed 
original and reduced systems. 

The data analyzed concerns 80 PE frames [7 of density p ~ 0.77— 0.78g/cm 3 . 
Each frame is generated by 8 PE chains of 1000 beads each, so finally our results 
on LK are based on 2240 distinct pairs of PE chains. 

First, we compute the normalized probability distribution of LK for pairs 
of PE chains in a frame and compare it to the LK for the corresponding re- 
duced pairs. Figure shows the probability distribution of the pairwise local 
periodic linking numbers for the original and CReTA reduced PE chains. The 
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LK 
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D^LK(open)-LK(closed) 




C 



LK 



D 



D^LK(cpen)-LK( closed) 



Figure 4: Figure A shows the normalized probability distribution of LK for the 
original chains before and after end-to-end closure and Figure B the probability 
distribution of the corresponding differences. Figure C shows the normalized 
probability distribution of LK for the reduced chains before and after end-to- 
end closure and Figure D shows the corresponding differences 



two distributions are quite similar. The mean absolute value of LK is 1.97 and 
the mean absolute value of LK of pairs of reduced chains is 2.01. Figure [3J3 
shows the probability distribution of the difference of LK before and after the 
application of the CReTA algorithm, i.e. D — LK or i g i na i — LK re d U ced- We note 
that the distribution is quite narrow (with standard deviation 0.30), has mean 
—0.19, and mean absolute difference (i.e. < \D\ >) of the local periodic linking 
between original and reduced chains of 0.29, indicating that LK is about the 
same for original and reduced pairs of chains. Thus, at least with respect to 
LK, one can study the entanglement of a PE melt in its reduced representation 
after the application of the CReTA algorithm. 

Next, we compute the local periodic linking number LK after performing 
end-to-end closure to the original chains and compare it to the LK after perform- 
ing end-to-end closure to the corresponding reduced chains. Figure[3P shows the 
normalized LK for the end-to-end closures of the original and the corresponding 
reduced chains. We notice a smaller overlap of the two distributions. The mean 
absolute value of LK for the original closed chains is 2.01 and the mean absolute 
value of LK for the reduced closed chains is 1.86. Figure [3p shows the normal- 
ized probability distribution of the difference D — LKorgCiosed — LKji e dciosed- 
We can see that the distribution is a bit broader having standard deviation 0.92 
in comparison with the open data. The mean absolute difference is 0.57. As we 
mentioned in Remark l3.41 we could expect LK to be invariant for closed chains 
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in most cases. The end-to-end closure is performed in this case before and after 
the application of the CReTA algorithm. Thus the observed differences between 
the original and reduced chains are most likely due to movements during the 
reduction algorithm in which a chain crossed the path of the end-to-end closure. 
We can see that the method of end-to-end closure retains much but not all the 
topological information from the original in the reduced systems. 

Next we compute the difference between the LK for open and closed orig- 
inal chains and the difference between the LK for open and closed reduced 
chains. Figures HR and HP show the normalized LK for the original open and 
the corresponding original closed chains and the normalized LK for the reduced 
open and the corresponding reduced closed chains respectively. In both cases, 
we notice an even smaller overlap between the two distributions. In Figures 
3J3 and HP is shown the normalized probability distribution of the difference 
LKo r gO P en - LK r g ciosed and of the difference LK Red0pe n - LK RedC i OS ed re- 
spectively. We can see that the distributions are broader, with standard devia- 
tions 2.00 and 1.91, respectively. The mean absolute difference is 1.57 and 1.52 
respectively. This shows that the local periodic linking number retains similar 
but distinct information for open and closed chains. 

5 Conclusions 

The local periodic linking number, LK, provides the means to effectively study 
the linking present in a periodic system as a measure of its entanglement. By 
employing LK and comparing the linking distributions before and after appli- 
cation of the CReTA algorithm we have shown that CReTA retains the linking 
information, and that LK reflects topological differences between open and 
closed chains. Our numerical results also show that the method of end-to-end 
closure of the chains is very sensitive to the movement of the chains, even if the 
TCs are preserved. This suggests that the local periodic linking number may 
be a stronger measure of entanglement for open chains. 
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